A stabilized linear finite element method for anisotropic poroelastodynamics with application to cardiac perfusion

We propose a variational multiscale method stabilization of linear finite element for nonlinear poroelasticity. Our approach is suitable the implicit time integration poroelastic formulations in which solid skeleton anisotropic and incompressible. A detailed numerical methodology presented monolithic formulation that includes both structural dynamics Darcy flow. implementation this verified using several hyperelastic benchmark cases, excellent agreement obtained with literature. Grid convergence studies hyperelastodynamics poroelastodynamics demonstrate second-order accurate. The capabilities our are demonstrated model left ventricle (LV) heart derived from human imaging data. Simulations indicate anisotropicity myocardium has substantial influence on pore pressure. Furthermore, temporal variations various components pressure (hydrostatic resulting changes volume fluid) correlated variation added mass LV, maximum being at peak systole. order magnitude good